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Abstract 

Under standard assumptions including stationary and serially uncorrelated Gaussian gravita¬ 
tional wave stochastic background signal and noise distributions, as well as homogenous detector 
sensitivities, the standard cross-correlation detection statistic is known to be optimal in the sense 
of minimizing the probability of a false dismissal at a fixed value of the probability of a false alarm. 
The focus of this paper is to analyze the comparative efficiency of this statistic, versus a simple 
alternative statistic obtained by cross-correlating the squared measurements, in situations that 
deviate from such standard assumptions. We find that differences in detector sensitivities have 
a large impact on the comparative efficiency of the cross-correlation detection statistic, which is 
dominated by the alternative statistic when these differences reach one order of magnitude. This 
effect holds even when both the signal and noise distributions are Gaussian. While the presence 
of non-Gaussian signals has no material impact for reasonable parameter values, the relative inef¬ 
ficiency of the cross-correlation statistic is less prominent for fat-tailed noise distributions but it 
is magnified in case noise distributions have skewness parameters of opposite signs. Our results 
suggest that introducing an alternative detection statistic can lead to noticeable sensitivity gains 
when noise distributions are possibly non-Gaussian and/or when detector sensitivities exhibit sub¬ 
stantial differences, a situation that is expected to hold in joint detections from Advanced LIGO 
and Advanced Virgo, in particular in the early phases of development of the detectors, or in joint 
detections from Advanced LIGO and Einstein Telescope. 


*Electronic address: lionel.martellini@edhec-risk.com 
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I. INTRODUCTION 


A stochastic background of gravitational-waves is expected to arise from the superposition 
of independent signals at different stages of the evolution of the Universe, that are too 
weak or too numerous to be resolved individually. This background can be of cosmological 
origin, from the amplihcation of vacuum fluctuations during inflation [2ll |22l 02] , pre Big 
Bang models [51051120], cosmic (super)strings |T2l [32l [361 01] or phase transitions [3^9], 
or of astrophysical origin, from sources since the beginning of stellar activity such as core 
collapses to neutron stars or black holes [SI 1221 ESI 0H], rotating neutron stars [SH ES] 
including magnetars [2S1 EH ESI 0S], phase transition [TB] or initial instabilities in young 
neutron stars PEI El 071 07] or compact binary mergers [21 [371 01 091 ED]. 

The detection of gravitational waves has become a question of central importance in as¬ 
trophysics and the detection of the cosmological contribution would have a profound impact 
on our understanding of the evolution of the Universe, as it represents a unique window on 
the very early stages up to a fraction of second after the Big Bang. An increasing range 
of efforts are dedicated to the design of improved detectors. The next generation of instru¬ 
ments (Advanced LIGO and Advanced Virgo [261 03]). which will start operating in 2015 
and 2016 respectively, are expected to be more than ten times more sensitive than their hrst 
generation counterparts. Besides, third-generation interferometers such as the European 
project Einstein Telescope (ET) [DD], currently under design study, are expected to fur¬ 
ther increase the likelihood of detecting the exceedingly small effects of gravitational waves. 
In parallel to the technological efforts towards the generation of sensitivity improvements 
for gravitational wave detectors, an increasing body of research is attempting to improve 
upon the efficiency of data analysis methodologies involved in stochastic gravitational wave 
background (SGWB) signals. 

The commonly used approach to the detection of GW stochastic background signals 
consists of cross-correlating the coherent measurements obtained from a pair of detectors. 
Under standard assumptions including stationary and serially uncorrelated Gaussian grav¬ 
itational wave stochastic background signal and noise distributions as well as homogenous 
detector sensitivities, the cross-correlation (GG) detection statistic is known to be optimal 
in the sense of minimizing the false dismissal probability at a hxed value of the false alarm 
probability (see for example [1102102]). Recent predictions based on population modeling 
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however suggest that, for many realistic astrophysical models, there may not be enough 
overlapping sources, resulting in the formation of a non-Gaussian background. It has also 
been shown that the background from cosmic strings could be dominated by a non-Gaussian 
contribution arising from the closest sources [I2l |36] . In the past decade a few methods have 
been proposed to search for a non-Gaussian stochastic background, including the probability 
horizon concept developed by m based on the temporal evolution of the loudest detected 
event on a single detector, the maximum likelihood statistic of [H] or [30] , which extends the 
standard analysis in the time domain in the case of parametric or non-parametric deviations 
of normality, the fourth-order correlation method from [30] , which uses fourth-order correla¬ 
tion between four detectors to measure the third and the fourth moments of the distribution 
of the GW signal, or the recent extension of the standard cross-correlation statistic by HU- 
While most of these papers maintain the assumption of Gaussian noise distributions so as 
to better focus on the impact of deviations from normality of the signal distribution, there 
is also ample evidence of strong deviations from the Gaussian assumption for noise distri¬ 
butions in gravitational waves detectors (see lasi), and relatively little is known about the 
impact of the presence of such non-Gaussian noise distributions on the efficiency of standard 
methods used for the detection of SGWB signals. Besides, the standard assumption that the 
two detectors have the same sensitivity may not hold strictly for joint observations by Ad¬ 
vanced LIGO and Advanced Virgo [26lll3], especially during the early stages of development 
of the detectors, or joint observations by Advanced LIGO and ET. 

The focus of this paper is to analyze the efficiency of the standard GG statistic in situa¬ 
tions that deviate from the aforementioned standard assumptions, and in particular involve 
deviations from the Gaussian assumption and/or the presence of detectors with heteroge¬ 
nous sensitivities. To do so we hrst introduce a simple alternative statistic obtained by 
cross-correlating the squared measurements, and we derive closed-form expressions for the 
mean and variance of this statistic as a function of the hrst four cumulants of the signal 
and noise distributions. We also show how to obtain consistent estimates for these parame¬ 
ters using a suitable extension of the likelihood function, for which we obtain an analytical 
expression. These results extend our previous results [30], where we have focussed on a 
situation involving a non-Gaussian signal distribution, but have maintained the assumption 
of a Gaussian noise distribution. Turning to a numerical analysis, we hnd that differences in 
detector sensitivities have a large impact on the comparative efficiency of the GG detection 
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statistic, which is dominated by the alternative statistic when these differences reach one 
order of magnitude. Remarkably, this effect holds even when both the signal and noise dis¬ 
tributions are Gaussian. While the presence of non-Gaussian signals has no material impact 
for reasonable parameter values, we hnd that the relative inefficiency of the GG statistic is 
less prominent in the presence of fat-tailed noise distributions, which imply an increase in 
the variance of the alternative detection statistic. On the other hand, the relative ineffi¬ 
ciency of the GG statistic is magnihed in case noise distributions have skewness parameters 
of opposite signs, a situation that leads to a reduction in the variance of the alternative de¬ 
tection statistic through a diversihcation effect. Overall, our results suggest that introducing 
an alternative detection statistic can potentially lead to noticeable sensitivity gains when 
noise distributions are non-Gaussian and/or when detector sensitivities exhibit substantial 
differences. 

The rest of the paper is organized as follows. In Section 2, we provide a brief review of 
the standard cross correlation statistic and introduce the alternative detection statistic. In 
Section 3, we perform a comparative analysis of the efficiency of the GG detection statistic 
versus the alternative detection statistic, and show that the latter dominates the former in 
a number of cases of potential practical relevance. In Section 4, we extend the maximum 
likelihood estimation techniques to a situation involving potentially non-Gaussian signal 
and non-Gaussian noise distributions so as to obtain consistent estimators not only for the 
variance but also the skewness and kurtosis of the signal and noise distributions, which are 
needed for implementing the alternative statistic. Finally, Section 5 contains a conclusion 
and suggestions for further research. 

II. INTRODUCING A NEW DETECTION STATISTIC 

In this Section, we hrst recall standard results related to the cross-correlation statistic. 
We then introduce an alternative statistic given by the cross-correlation of squared measure¬ 
ments. 
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A. Assumptions and Notation 


Consider two gravitational wave detectors. The output of each detector is a collection 
of dimensionless strain measurements. Suppose that N such measurements are made by 
each detector at regular time intervals. Denote these measurements by a T x 2 matrix h 
with components h^, where i = 1,2 labels the detector, and t = 1,2, ...,iV is the discrete 
date of measurement. To determine whether or not the data h contains some desired signal, 
one usually compares the value of some detection statistic T (h) to some threshold value T*. 
If T (h) is greater than the threshold value T*, one concludes that a signal is present and 
otherwise one concludes that no signal is present. A detection statistic is said to be optimal 
if it yields the smallest probability of mistakenly concluding a signal is present (probability 
of a false alarm, or pfa) after choosing a threshold which hxes the probability for mistakenly 
concluding that a signal is absent (probability of a false dismissal, or pfd). 

We hrst decompose the measurement output for detector i in terms of noise versus signal, 
which gives when written in terms of random variables: 

'Hi= Mi+Si 

where A/i denotes the noise detected by the detector i and Si denotes the signal detected 
by the detector i so that Hi is the total measurement for the detector i. If we now assume 
that the detectors are coincident and coaligned (i.e., they have identical location and arm 
orientations), we obtain that the signal received by both detectors is drawn from the same 
distribution. Under this assumption, we would have that: 

Si = S2 = S 

In terms of the realization of such random variables for either one of the two detectors, 
we note: 

T Sf 

Given that both signal and noise distributions can potentially be non-Gaussian, we de¬ 
noted by Cj, j = 1,2,3,4, the hrst four cumulants of the signal distribution, and by qj, 
j = 1,2,3,4, the hrst four cumulants of the noise distribution for detector i, with i = 1 or 

2 . 
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Let us recall that for a random variable X with density function denoted by fx (here 
X = 5 or A/i, for i = 1, 2), we can introduce the moment generating function: 


Mx (t) = E [e*^] = / (a:) dx 


( 1 ) 


which is related to the characteristic distribution ipx, i-e., the Fourier transform of the 
function fx, by ipx {t) = Mx {it). The (non central) moment of the distribution of 
the random variable X is given by the derivative of the moment-generating function Mx 


taken at t = 0 (hence the name moment generating function): fij = M^^ (0) = {—i)^ fjx’ (0). 

OO 

Using the Taylor expansion of the exponential function around 0, e* = V —, we obtain a 

j=o r- 

new expression for the characteristic function: 




J ^—n J * 


( 2 ) 


j=o J=0 

We also introduce the cumulant generating function mx as the logarithm of the moment 
generating function: 

“ (t)^ 


mx {t) = log Mx (t) = log ^ 




J'- 


(3) 


A Taylor expansion of the cumulant generating function mx is given by a series of the 
following form: 


P 


mx {t) = mx (0) ^ (0) 

i=i 


(4) 


and we dehne cj = m^x (0) the cumulant of the random variable X.A moments-to- 
cumulants relationship can be obtained by expanding the exponential and equating coeffi¬ 
cients of P in: 


Mx (t) = exp [mx (t)] 


^ P 

2^ = exp 

j=o 


P 




j=i 


J'- 




(5) 


Conversely, a cumulants-to-moments relationship is obtained by expanding the logarith¬ 
mic and equating coefficients of P in mx {t) = logMx (t). In particular we have: 


Cl = m{ (0) = Pi = /i (6) 

C 2 = m" (0) = P 2 - /ii = 0 -^ (7) 

C3 = mf^ (0) = P3 - 3p2yUi + 2p? (8) 

C 4 = (0) = p4 — 4p3pi — 3p2 + 12p2hi “ 6 pi (9) 
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We note that the first cumulant is equal to the first moment (the mean), and the second 
cumulant is equal to the second-centered moment (the variance). For the Gaussian distri¬ 
bution with mean /i and variance we have Ci = fi, C 2 = cr^, and = 0 for k > 2. 
This allows us to identify deviations from the Gaussian assumption through the presence 
of non-zero 3rd- and 4th-order cumulants, C 3 and C 4 , which are sometimes normalized so as 

^3 

to transform into skewness and kurtosis parameters, respectively dehned as: skw = - 1 ^ 

^2 

and kurt = ^. In our application, it should be noted that signal and noise distributions 
are centered and therefore we have Ci = Ci^i = C 2 ,i = 0. We also use the notation C 2 = 

Cl ,2 = and C 2,2 = where a, ai, and (T 2 denote the standard-deviations for the signal, 

detector 1 and detector 2 distributions, respectively. 

B. Distribution of the Cross-Correlation Detection Statistic 

We hrst define the standard cross correlation detection statistic DScc as: 

1 ^ 

DS^ = -J2^itn2t ( 10 ) 

t=i 

where 'Hu = Mu -|- 5*, for i = 1 , 2 , and where Mu and for 1 < t < T, are T independent 
copies of the random variables Mi and S, respectively. 

We have: 

/ T T T ^ \ 

DScc = 7^ I A/ifA /21 + '^^MitSt + A/2t>Si -|- I (11) 

\ 1=1 1=1 1=1 1=1 / 

A signal is presumed to be detected when the detection statistic DScc exceeds a given 

detection threshold DT : 

DScc > DT (12) 

We typically select the detection threshold DT such that pfa = x%, for a given conhdence 
level x%, where the probability of a false alarm is given by the probability to exceed the 
threshold in a situation where there is no signal: 

pfa = PriDScc > DT\Hu = Mu) (13) 

Obviously, pfa is independent of the signal distribution. What depends on the signal 
distribution is the probability of a false dismissal pfd given by the probability that the 



detection statistic remains below the threshold even if there is a signal: 


pfd = FriDS,, < DT\ Uu = Mu + St) 


( 14 ) 


By the central limit theorem, it can be shown that the detection statistic DScc is asymp¬ 
totically normally distribnted whether or not the signal and noise distributions are Gaussian 


(see [20] for more details in the case of a non-Gaussian signal). In this situation, the distri¬ 


bution of the GG detection statistic is fully characterized by its mean and variance, which 
can be explicitly obtained as follows: 



-T [Yar {M 1 M 2 ) + Yar (MiS) + Yar {M 2 S) + Yar (5^)] 

(Yar (Ml) Yar {M 2 ) + Yar (Mi) Yar (5) + Yar {M 2 ) Yar {S) + E (5^) - (E 
(cr^cTg + afM -|- + C 4 + SM — = — [cfal + a\c? Y (y\c? Y ‘iM Y C 4 ) 


Hence, we obtain that for general signal and noise distributions the cross-correlation 
detection statistic DScc is asymptotically normally distributed, with mean o? and variance 



It should be noted that the variance of the 


detection statistic is identical whether or not the noise distributions are Gaussian since it 
does not depend on the higher order cumulants of the noise distributions, while it depends 
on the fourth-order cumulant of the signal distribution. This is because the standard cross¬ 
correlation detection statistic involves the squared value of the signal distributions, while 
the noise distributions are not squared. In what follows, we discuss the introduction of a 
new detection statistic that would make the detection procedure explicitly dependent upon 
the higher order cumulants of the noise distribution, and which will be found to dominate 
the cross correlation statistic for some realistic parameter values. 

C. Introducing an Alternative Detection Statistic 

For a Gaussian signal, the cross-correlation detection statistic can be shown to be optimal 
in the sense of minimizing the false dismissal probability at a hxed value of the false alarm 
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probability, a result which holds under restrictive assumptions [H] including stationary and 

serially uncorrelated Gaussian gravitational wave stochastic background signal and noise 

distributions. In the general non-Gaussian case, the cross-correlation detection statistic 

may not be optimal, and may be dominated by an alternative detection statistic, which can 

1 ^ 

be written in general as some function of the observations / (Hii, ^- 21 ) 1=1 t '^ 11 ^ 21 - 

It is unclear how one could derive an optimal detection statistic in a fully general setting, 
and we introduce in what follows a simple heuristic alternative detection statistic, denoted 
by DSait, which is given by the cross-correlation of squared detector measurements: 



(15) 


t=l 


By the central limit theorem, we know again that DSait is asymptotically Gaussian, and 
the asymptotic distribution for the detection statistic is therefore fully characterized by its 
hrst two moments Halt and We hrst have: 



(16) 


(17) 


22, 22, 22io4i 

(7^(72+ (TiOl +(72(7 + oa + C 4 


(18) 


In case the signal is absent (c 4 = = 0), the expression for the mean value of the 


alternative detection statistic (denoted by Halt ^^is case, where ns stands for no signal) 
further simplihes into: 



(19) 


.2 

alt 




We then compute a[ 




t=i 



t=i 


^Var [(Ari + 5)'(A/'2 + 5f] 


We note that: 


Var [(M + 5)' (W 2 + 5)"] = E [(M + 5)" (A4 + 5)^] - (E [(M + 5)' (A4 + 5)"])' 
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After some tedious computations, we first obtain: 

E [(A/l + S) (A /2 + iS) ] = fj,s + 6 /r 6 /^l ,2 + 6 /i 6 /i 2,2 + d/Xs/ii^S + 4p,5/i2,3 + At 4 Ail ,4 + 6 /i 4 / 42,4 

+36/i4/4i,2/^2,2 + 24/43/ii^3/i2,2 + 24/i3/4i^2/^2,3 + 6/i2/^l,4/^2,2 + 6/42/^!,2/^2,4 
+ 16/42/^1,3/^2,3 + /^l,4/^2,4 (20) 

where we use the following notation for the higher order moments of the signal and noise 
distributions for j = 1, 

= E [S^] 

Hij = E [A/7] 

H2,j = E [A/7] 

Finally we obtain: 

Tail, = E [(M + ‘5)' (A4 + 5)"] - (E [(A^i + 5)' (A4 + 5)'])' 

which after more tedious calculation gives the following expression for the variance of the 
alternative detection statistic: 

^alt — y (A'-S — /^4 + 2 (3/i6 — /i2h4) (hi,2 + ^ 2 , 2 ) + 4/i5 (hi,3 + h2,3) 

+h4hl,4 ~ h2hl,2 + h4h2,4 ~ h2h2,2 + 2 (l7/i4 — /i^) hl,2h2,2 
+24/i3 (hl,3h2,2 + hl,2h2,3) + 16h2hl,3h2,3 + 2/i2h2,2 (3hl,4 “ hi,2) 

+2h2hl,2 (3h2,4 — h2,2) + hl,4h2,4 “ hl,2h2,2) (21) 

If we now assume a symmetric signal distribution (/ii,/i 3 ,/is,/i 7 = 0), the expression for 
the variance of the detection statistic simplihes into: 

X 

^llt = y (hs — h4 + 2 (3h6 — h2h4) (hl,2 + h2,2) (22) 

+h4hl,4 ~ h2hl,2 + h4h2,4 ~ h2h2,2 + 2 (l7/i4 — /i^) hl,2h2,2 
+ 16/i2hl,3h2,3 + 2/i2h2,2 (3hl,4 ” hi, 2 ) 

+2h2hl,2 (3h2,4 ~ h2,2) + hl,4h2,4 ~ hl,2h2,2) (23) 

If the noise distributions are also symmetric (hi, 3 )h 2,3 = 0), the expression further sim- 
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plifies into: 


^llt = y (A'-S - /^4 + 2 (3/i6 - ^^2|^4:) (/il,2 + /^2,2) 

+/i4/il^4 — 2 + Ai4Ai2,4 ~ fAl^2,2 + 2 [l7f^4 — 1^1,21^2,2 

+2/i2/^2,2 (3/il,4 ~ /^1,2) + 2/i2/^l,2 (3/^2,4 ~ /^2,2) + /^l,4/^2,4 ~ 2/^2,2) (24) 

In case the signal is absent (p-g = /ie = h 4 = /^2 = 0), the expression for the variance of 
the alternative detection statistic (denoted by becomes: 

f (/^I,4/42,4 - /4?, 2 / 4 ^ 2 ) = ^ ((^1,4 + 3a^) (c2,4 + 3a^) - afcT^) (25) 

Clearly, the variance of the alternative detection statistic is higher when the higher order 
cumulants of the noise are not zero, which will have implications for the sensitivity of the 
detection procedure. The higher variance of the distribution of the alternative detection 
statistic in the non-Gaussian case (that is the case when noise distribution are potentially 
non-Gaussian) implies that the alternative detection statistic has fatter tails compared the 
Gaussian case (that is the case when noise distributions are Gaussian). As a result, the 
detection threshold corresponding to a given pfd will be lower in the non-Gaussian case, 
which in turn allows for the detection of fainter signals when the presence of non-Gaussianity 
is taken in to account with respect to a situation where the observer uses the alternative 
detection statistic while wrongly assuming that the underlying signal and noise distributions 
are Gaussian. In other words, for an observer using the alternative detection statistic, 
taking into account the non-Gaussianity of the signal and noise distributions will improve 
the detection methodology. 

An outstanding question, however, remains with respect to whether or not the observer 
would be better off using the alternative versus the standard cross-correlation statistic. 
While the alternative detection statistic has no claim to optimality, and is expected to be 
dominated by the cross-correlation statistic when signal and noise distributions are Gaussian, 
it may in principle dominate the standard cross-correlation statistic in the non-Gaussian case 
since the optimality of the GG statistic has not been established in this more general setting. 
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III. COMPARATIVE EFFICIENCY OF THE CROSS-CORRELATION VERSUS 
THE ALTERNATIVE DETECTION STATISTIC 


To compare the performance of the standard and alternative detection statistic in the 
general non-Ganssian case, we nse the following mnlti-step procednre. 


• Step 1: We select a set of parameter valnes for the signal and noise distribntions, and 
we apply the transformation that allows one to tnrn the cnmnlants into corresponding 
moments, which will be needed in the expressions for the mean and variance of the 
alternative detection statistic. Assnming for simplicity symmetric signal and noise 
distribntions, we have: 

/i2 = C2 = Q;^,/ii_2 = Ci^2 = = C2,2 = <^2 (26) 

/i4 = C 4 + 3C2, /il,4 = Ci^4 + 3c^ 25 h2,4 = C2,4 + 3C2 2 (27) 

he ~ Cg + I 5 C 4 C 2 + 15c2 (28) 

/is = C8 + 28 c6C2 + 35 C 4 + 2 IOC 4 C 2 + 105c2 (29) 


We may then obtain the corresponding valnes for fiait and in the presence of the 


signal, nsing equations (16) and (24), as well as the corresponding valnes for and 


in the absence of the signal, nsing eqnations (19) and (25). 


• Step 2: We select a given pfa valne, taken in the nnmerical analysis that follows to 
be 5%, and obtain the corresponding thresholds for both the standard and alternative 
detection statistics, denoted respectively by DScc and DSait (which are fnnctions of 
the selected pfa valne), nsing the following eqnations: 

pfa = Pt{DS,,> DT,, {pfa) \ = Mt) (30) 

= PiiDSait > DTait {pfa)\Uu = Mt) (31) 


based npon the following Ganssian distribntions for the detection statistic in case the 
signal is absent: 


T 

t ^ V (^0, ^a'ia'i) 


t=l 

T 


DTait = ^ ^ ^ ((ci,4 + 3a^) (c2,4 + 3a^) - 

t=i ^ 


where we take T = 10^ in the base case 
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• Step 3: We compute the probability of a false dismissal corresponding to the standard 
cross-correlation statistic and the probability of a false dismissal corresponding to the 
alternative detection statistic using: 

pfdcc = Pr (DS'cc < DT^c = Uit + St) (32) 

pfdait = 'PT{DSait < DTait ipfa)\'Hit = Mit + St) (33) 

based upon the following Gaussian distributions for the detection statistic in case the 
signal is present: 

DScc = ^ 'Hit'H2t A/" fa"^, ^ + 2a^ + C4)) 

t=i ^ ' 

1 ^ 

DT^it = 

t=l 

If we can find a set of parameter values and a pfa value such that pfdait < pfdcc, we would 
then prove that the standard cross-correlation analysis is not always optimal, and we would 
also show that the alternative statistic we have introduced allows for a more efficient detec¬ 
tion procedure, at least for the selected set of parameter values. In what follows, we analyze 
the probability of a false alarm for both statistics in various situations involving homoge¬ 
nous versus heterogenous detector sensitivities, as well as various assumptions regarding 
the higher order moments of the noise distributions. Note that the signal distribution is as¬ 
sumed to be Gaussian in all the result that we present below. In unreported results, we have 
analyzed the relative efficiency of the two statistics in situations involving a non-Gaussian 
signal, and have found only very small differences with respect to the Gaussian signal case. 
Indeed, the signal is assumed to be small compared to the noise in realistic situations, and 
therefore the impact of deviations from the Gaussian assumption at the signal level will be 
dwarfed by the impact of deviations from the Gaussian assumption at the noise level. 

In Fig. we hrst plot the probability of a false dismissal (pfd) as a function of the 4th 
cumulant of the noise distribution, assumed to be identical for both detectors (ci ^4 = 02 , 4 ), 
over a resonnable range of values expected for classical distributions such as Gaussian, 
Hypersecant, Logistic or Laplace (see 1301 ), for the GG statistic (in blue) and the alternative 
statistic (in red and green). Here we assume that the signal is Gaussian (cs = C 4 = 0) 
and that the noise distributions are symmetric (ci^a = 02,3 = 0). The parameter a is 
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TABLE I: Noise variance levels ai and cj 2 , and ratio ri 2 = —calculated as dfSn(f), 

where /min = 10 Hz and /max = 250 Hz is the typical frequency band used for the cross-correlation 
analysis, for Advanced LIGO with Advanced Virgo (aLIGO and AdV) at design sensitivities and 
during the early phases of development of the detectors (early, middle and late), and for Einstein 
Telescope (ET-D sensitivity) with LIGO Red [16], a possible Advanced LIGO sensitivity upgrade. 


Pair 


^2 

n2 

AdV - aLIGO 

3.6 X 10“^“^ 

1.7 X 10“^2 


early- middle (6 months) 

3.8 X 10“^^ 

1.0 X 10"‘^3 

371 

middle - late (9 months) 

1.5 X 10"^^ 

3.6 X 10"^^ 

402 

late - design (12 months) 

3.2 X 10"^2 

3.6 X 

88 

LIGO Red - ET-D 

6.0 X 10"^^ 

1.9 X 

31 


chosen so that the signal-to-noise ratio is SNR = \/T - = 3.28, a value yielding a 

cri(J2 

probability of false alarm pfd = 5% for the CC statistic in the homogeneous case where 
detectors 1 and 2 have the same sensitivity, i.e. when the ratio between the detector noise 

( 7 ^ 

variances is ri 2 = ^ = 1. The various red plots indicate different values of this ratio, ri 2 = 
1,10,100, where detectors 1 and 2 are chosen so that ri 2 > 1. The green plots correspond 
to realistic values for the cross-correlation between Advanced LIGO and Advanced Virgo 
at their nominal sensitivity (ri 2 = 48), and between Einstein Telescope and LIGO Red 
[To] , a possible Advanced LIGO sensitivity upgrade (ri 2 = 30). We also considered a value 
ri 2 = 400 corresponding to the maximum expected cross-correlation between Advanced 
LIGO and Advanced Virgo during the early phases of the development of the detectors [T]. 
The projected nominal and early sensitivities, in term of the square root of the power spectral 
density Sn, of Advanced LIGO and Advanced Virgo [26],|l3], along with the LIGO Red noise 
curve [16] and the proposed Einstein Telescope sensitivity ET-D [23] are plotted on Fig. 
The corresponding noise variances, calculated as dfSn(f), where /min = 10 Hz 

and /max = 250 Hz is the typical frequency band used for the cross-correlation analysis [31] 
are reported in Table 

We confirm that pfdcc = 5% when the two detectors have equal sensitivity (ri 2 = 1) and 
we hnd that pfdcc < vfdait in this case, as expected. On the other hand, as we let the ratio 
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ri 2 increase, we find that the probability of a false dismissal for the CC statistic decreases very 
slightly, while the probability of a false dismissal for the alternative statistic decreases very 
fast. As a result, we have that pfdcc > pfdait when ri 2 >> 10. Note that pfdait is almost 
0 when ri 2 = 100, a situation in which the signal appears large compared to the noise in the 
most sensitive detector. We therefore obtain that the standard CC statistic can be dominated 
by the alternative statistic when detector sensitivities exhibit substantial differences, even 
when both signal and noise distributions are Gaussian. In fact, the domination of the 
alternative statistic decreases as Ci ^4 = C 2,4 increases. This can be explained by the fact that 
increases in ci ^4 and C 2,4 lead to increase in the variance of the alternative detection statistic, 
which is detrimental to the performance of the detection methodology. We notice that for 
the values of ri 2 expected for the current and next generations of detectors, the alternative 
statistic usually performs better than the CC statistic, except for the most pessimistic case 
when ci ^4 = C 2,4 > 8 for the cross-correlation between LIGO Red and ET (ri 2 ~ 48). 

In Fig. 1^ we consider the dual perspective, where the probability of a false dismissal is 
plotted against the ratio ri 2 for typical values of Ci ^4 = C 2,4 = 0, 3 and 10. Again, as the 
ratio increases, we hnd a small deterioration in the performance of the GG statistic (the 
probability of false dismissal increases from pfdcc — 5% for ri 2 = 1 to pfdcc — 8% for 
ri 2 = 1000) and a substantial improvement for the alternative statistic. In the case of a 
normally distributed noise process (ci ^4 = C 2,4 = 0), the alternative statistic outperforms the 
GG statistic for values ri 2 > 6 and the probability of a false dismissal becomes negligible, 
as small as pfdait ~ 4 x 10“^° for the pair Advanced LIGO and Advanced Virgo (ri 2 — 48) 
and pfdait ~ 1 x 10“® for the pair LIGO Red and ET (ri 2 — 30), which translates into 
a gain of 8 and 4 orders of magnitude, respectively, compared to the GG statistic. For 
the cross-correlation between Advanced LIGO and Advanced Virgo during the early stages 
of development (ri 2 — 400), the probability of a false dismissal is almost zero {pfdait ~ 
2 X lO-^'^^)^ 

So as to better understand why the presence of heterogenous detectors has such a strong 
impact on the relative efficiency of the GG versus alternative detection statistic, we consider 
two contrasted situations, a homogenous case situation (Gl), with cxi = (72 = a and a 
heterogenous case situation (G2), with ai = — and (72 = We note that by construction 

(7 

the product (7^(7| = in both cases but the sum + is different and substantially higher 
in the heterogenous case situation. 
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FIG. 1: Comparative efficiency of the CC and alternative detection statistics. We take T = 10® 

and pfa = 5% and we assume that the signal is Ganssian (ca = C 4 = 0) and that the noise 

distributions are symmetric (ci^a = 02,3 = 0). Various red plots correspond to different values for 
2 

the ratio ri 2 = ^ = 1,10,100,1000. The parameter a is chosen so that the signal to noise ratio 

r- 

SNR = vT- = 3.28, a value that yields a pfd = 5% for the CC statistic in the homogeneous 

< 72(72 

case ri 2 = 1. When ci ,4 = C 2,4 = 0, pfdait — 0.2,2 x 10“^°, 3 x 10“^° for ri 2 = 1,10,100. The 
green plots correspond to realistic values for the cross-correlation between Advanced LIGO and 
Advanced Virgo at their nominal sensitivities (ri 2 = 48) and between Einstein Telescope and 
Advanced LIGO with the LIGO Red possible upgraded sensitivity (ri 2 = 30). We also considered 
an average value of ri 2 = 400 corresponding to the cross-correlation between Advanced LIGO and 
Advanced Virgo during the early phases of the development of the detectors. When 04^4 = 02,4 = 0, 
pfdait - 1 X 10“®, 4 X 10-1°, 3 X 10-20,2 x 10-^®® for ri 2 = 30, 48 and 400. 
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f(Hz) 


FIG. 2: Expected sensitivity of Advanced LIGO and Advanced Virgo (blue and black continuous 
lines), LIGORed (dashed red line) and ET-D (green continuous line). The evolution of the sensi¬ 
tivity during the Advanced LIGO and Virgo early, middle and late phases are also shown in dashed 
blue and black lines. 


Assuming Gaussian distributions for signal and noise, we have: 




= 0, yielding the same value in Cl and C2 


Raft 


2,ns 


= crfa 2 , yielding the same value in Cl and C2 




^cc - 


2,ns _ 

^alt 


-, yielding the same value in Cl and C2 




, yielding the same value in Cl and C2 
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FIG. 3: Probability of a false alarm for the CC and alternative detection statistics as a function of 
the ratio of detector sensitivities. We take T = 10^ and pfa = 5% and we assume that the signal 
is Gaussian (cs = C 4 = 0) and that the noise distributions are symmetric (ci^s = 02,3 = 0). Various 
red plots correspond to different choices for the pair of parameters C 14 and 02 , 4 . The parameter a 
is chosen so that the signal to noise ratio SNR = \/T-= 3.28, a value that yields a pfd = 5% 

0'2(T2 

2 

for the CC statistic in the homogeneous case ri 2 = —| = 1. 


and we also have: 


Rcc 

f^alt 

^2 

^cc 

2 

^ alt 


yielding the same valne in Cl and C 2 
cr^(T 2 + (di + CT 2 ) + + C 4 , yielding a greater valne for C2 

— {dial + {dl + dl) + 2a'^ + C 4 ) , yielding a greater valne for C2 


expression given in Eq. 24, yielding a greater valne for C2 
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As a result, we find that the probability of a false alarm is the same in Cl and C2 for 
both the standard CC and the alternative detection statistics, since this probability only 
depends on the distributions of the detection statistics in the absence of a signal, which are 
the same for Cl and C2. In the presence of a signal, we find for the alternative statistic that 
moving from the homogenous case (Cl) to the heterogeneous case (C2) leads to an increase 
in the mean, which has a positive impact on sensitivity, and an increase in the variance, 
which has a negative impact on sensitivity. Overall, the net effect is positive, as can be seen 
from Fig. |^and|^ For the cross-correlation statistic, the mean is not impacted but there is 
an increase in variance, which is detrimental to detection. Overall, we conhrm that the case 
with heterogeneous sensitivities is more favorable for the alternative statistic than it is for 
the CC detection statistic. 

We now turn to the analysis of the impact of the 3rd moment of the noise distribution, 

which has been assumed to be zero so far. In Fig. we show the probability of a false 

dismissal for the CC statistic (in blue) and the alternative statistic (in red) as a function 

of the third-order cumulant of the distribution of the noise for the hrst detector. We take 

T = 10® and pfa = 5% and we assume that the signal is Gaussian (cs = C4 = 0). We consider 

(^2 

the homogeneous case ri2 = ^ = 1. As usual, the parameter a is chosen so that the signal- 

/— Q!^ 

to-noise ratio SNR = yT -= 3.28, a value that yields a pfd = 5% for the CC statistic 

0 ' 20'2 

for ri2 = 1. Various red plots correspond to different values for 02,3, with C4,j = 0 for i = 1, 2, 
except for the crossed red line, where it is equal to 1. We hnd that the alternative statistic 
may dominate the CC statistic when ci^s and 02,3 are of opposite signs. Indeed, the 3rd 
higher-order cumulants of the noise distributions do not impact the mean and variance of the 
CC statistic, but have an impact on the variance of the alternative statistic. Unlike the 4th 
order cumulant that is always positive, the 3rd higher-order cumulants can principle take on 
any negative or positive value, and taking them of opposite signs leads to a reduction in the 
variance of the alternative detection statistic. Indeed, these parameters enter the expression 
for the variance of the alternative statistic in Eq. through the term -l-16/r2/ii,3yU2,35 which 
is negative when noise distributions have skewness parameters of opposite signs, suggesting 
a diversihcation effect (note that ps = C3 for centered distributions). On the other hand, we 
note again that an increase in (crossed red line in Fig. leads to a deterioration of the 
performance of the alternative statistic. 
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FIG. 4; Impact of the skewness of the detector noise on the relative efficiency of the CC and 

alternative detection statistics for homogeneous detector sensitivities.. We take T = 10^ and 

pfa = 5% and we assume that the signal is Gaussian (ca = C 4 = 0). We consider the homogeneous 
2 2 
case ri 2 = —= 1. The parameter a is chosen so that the signal to noise ratio SNR = y/T -= 

(T 2 (72 0-2 

^2 

3.28, a value that yields a pfd = 5% for the GC statistic for ^ = 1 . Various red plots correspond 
to different choices of value for ci ,3 and 02 , 3 , with C 4 ,j = 0 for i = 1 , 2 , except for the crossed red 
line, where it is equal to 1 . 


In Fig.[^ we repeat the analysis but focus on the case of heterogenous detector sensitivities 


by taking ri2 = ^ = 10, while we had ri2 = 1 in Fig. 




We hnd again that the alternative statistic domina 


es the CC statistic, an effect that is 


stronger when ci,3 and 02,3 are of opposite signs, but decreases when C4,i increases. 
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FIG. 5; Impact of the skewness of the detector noise on the relative efficiency of the CC and 

alternative detection statistics for heterogeneous detector sensitivities. We take T = 10^ and 

pfa = 5% and we assume that the signal is Gaussian (ca = C 4 = 0). We consider the homogeneous 
2 2 

case ri 2 = —^ = 10. The parameter a is chosen so that the signal to noise ratio SNR = y/T -= 

a 2 U2cr2 

3.28, a value that yields a pfd = 5% for the CC statistic for ri 2 = 1. Various red plots correspond 
to different choices of value for ci^a and 02 , 3 , with 04 ^* = 0 for i = 1 , 2 , except for the crossed red 
line, where it is equal to 1 . 

IV. ESTIMATION METHODS FOR NON-GAUSSIAN SIGNAL AND NON- 
GAUSSIAN NOISE DISTRIBUTIONS 

The analysis in the previous Section suggests that the use of the alternative statistic may 
lead to noticeable sensitivity gains when noise distributions are non-Gaussian and/or when 
detector sensitivities exhibit substantial differences. It should be noted, however, that using 
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the alternative statistic requires in the non-Gaussian case the use of robust estimates not 
only for the variance, but also the skewness and kurtosis, of the signal and noise distributions. 
In what follows, we show how to obtain such estimates by extending standard maximum 
likelihood estimation methodologies to situations involving possibly non-Gaussian signal and 
noise distributions. As such, these results generalize early results reported by 133 who have 
focussed on a situation involving a non-Gaussian signal distribution, while maintaining the 
assumption of a Gaussian noise distribution. 

We denote by /„. and fs, respectively, the density function for the noise and the signal: 

Pr {Afi e [n,n + dn]) = (n) dn 

Pr {S G [s, s -|- sn]) = fs (s) ds 


We also denote by /„ = /„ {nu, '^2t)f=i t joint probability distribution for the noise 
in the two detectors. The standard Bayesian approach for signal detection consists in hnding 
the value for the unknown parameters so as to minimize the false dismissal probability at 
a hxed value of the false alarm probability. This criteria, known as the Neyman-Pearson 
criteria, is uniquely dehned in terms of the so-called likelihood ratio A given by: 


A 


Ph\x=l 

Ph\x =0 


where Ph\x=i (respectively, Ph\x=o) conditional density for the measurement output 

if a signal is present (respectively, absent). A natural approximation of the likelihood ratio 
is the maximum likelihood detection statistic dehned by 


Amt. — 


max / fs\x=i (s) fn\x=i (^ “ s) ds 

,C72 


max /„! V.0 W 


(34) 


<Ti,cr2 


and the maximum likelihood estimators for the unknown signal and noise standard-deviation 
parameters a, ai and a 2 are given as the corresponding likelihood maximizing quantities. 


Full Gaussian Case 


It is typically assumed that both the noise and signal are normally distributed, /„ and 
fs are Gaussian probability distribution functions, that is we may assume: 


frii iP-it) 




\/^i 


for i = 1,2 


TTcr,: 


1 (h-pP 
fs{st) = ,— e 2-" 




7ra 
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where we also assume that both the noise and signal are weakly stationary processes so that 
their moments are constant through time. We further assume the noise in detector one and 
two are uncorrelated with a zero mean for both detectors. Under these assumptions, we 
have: 


fn {nit,n2t) 



2 naia2 


and hnally, assuming zero serial correlation: 


fn — fn (^It) 


1 


n-i . rin 


t=i 


27r(Ticr2 


^ .^ “ n 


(H-nr 
-e 2 a'^ 


Typically, one also assume that the mean value (3 for the signal is zero so that the unknown 


parameters are a, cxi and <72. Then, we have that the denominator of equation (34) is given 
by: 

T 


max fn\x=o {h) = maxTT 






27rCTiCT2 


'“It ^2t 


-e ^"^1 ^‘"2 = max- 


0-1 ,<5-2 (27rCTiCT2)" 


exp 


Introducing for z = 1, 2: 


^ 2cr2 ^ 2 a^ 
t=l 1 i=l 2 


a? = ^ li: 


T 


t=i 


we hnally have that: 


max/„ = max-? 

cri,(T2 0-1,(T2 (27rcri(T2)" 


exp 


2 


(35) 


It is straightforward to see that the maximum for equation (35) is reached for af = a‘1 
and that maximum is given by: 


max/„ = -— _ j, exp 

0-1,0-2 (27r(JiCr2) 


T 

■ 2(1 + 1 ) 




Finally, we obtain the following expression: 


Aatl — 


max / fs (s) fn {h - s) ds 

a,(Ji ,(J2 

max/„ 

^1,0-2 


^ P + OO 

= (27rcfia2)^ exp (T) max IT / 


fs (St) 


27iaia2 


exp 


{hit - Stf {h2t - Stf 


2a? 


2 a? 


dst 


= max 


n 


7;; ^ /■+00 

<7ia2 / 


>1^1’1^2 *^1^2 J -00 


fs (st) exp 


(hit - stf {h2t - stf 


2 af 


2 al 


+ 1 


dst 
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We now specialize the analysis of the specihc situation where the signal has a Gaussian 
distribution. In this case, and maintaining the assumption that the mean value for the signal 
is zero, we obtain: 


/. - /. - n 


-e 


We thus have: 


^ML = max TT 

Q,(Tl,C72 


cri(J2 


^+oo 


\/‘2,710.(71(72 J-c 


exp 


Si {hit - {h2t - Stf 


2a2 


2 al 


2 al 


+ 1 


dst 


After tedious calculations, one obtains (see my- 


Aml = max 


7172 


0-1,0-2>0 j ^yofoY+~7fa^~+~o^a^ 


exp 


I ^2 I 2q;2 
^ 1^2 


2 ( ^ + ;4 + ^ I ^7 2a: 


^2 


(36) 


One can show that the maximum is reached for: 


2 ^2 _ /—2\ + 

o = a = [a ) 

2 ^2 _ /—2 ^2\ + 

= [^i - ) 

where (x)^ = a: if x > 0 and (x)"'" = 0 otherwise, which arise because of the positivity 
constraints on a, 7i and <72 in he maximization procedure. 

The corresponding detection statistic is: 

/ ^4 \ —T/2 

= 1 - ^2 (37) 

\ ^ 1 ^ 2 / 

The cross-correlation statistic can be obtained from A^^ via a monotonic transfor¬ 
mation which preserves false dismissal versus false alarm curves (see again n): 


B. Gaussian Signal and Non-Gaussian Noise 

In [30] , the Gaussian assumption was maintained for the detector noise distribution, but 
relaxed for the signal distribution. In what follows, we consider the opposite situation, 
namely a normally distributed signal, and a potentially non-Gaussian noise distribution. In 
other words, we assume: 
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1 


fs{St) = 

frii {nu) 7 ^ 


y/^a 

1 


e ^ 


2o- 


i for i = 1, 2 


As in [30] , we propose to use a semi-parametric approach which allows one to approximate 
the unknown density as a transformation of a reference function (typically the Gaussian 
density), involving higher-order moments/cumulants of the unknown distribution. This 
approach has been heavily used in statistical problems involving a mild departure from the 
Gaussian distribution. In what follows, we will show that it allows us to obtain an analytical 
derivation of the nearly optimal maximum likelihood detection statistics for non-Gaussian 
gravitational wave stochastic backgrounds. 

We want to approximate fm, the density function of the unknown distribution of the noise 
distribution A/), as a function of the Gaussian density function f^. (x) and a multiplicative 
deviation gn^ (x) from the Gaussian density function. To achieve this objective, we use the 
Edgeworth expansion^ which is based on the assumption that the unknown signal distribution 
is the sum of normalized i.i.d. (non necessarily Gaussian) variables. In other words, it 
provides asymptotic correction terms to the Gentral Limit Theorem up to an order that 
depends on the number of moments available. When taken to the fourth-order level, the 
Edgeworth expansion reads as follows (see for example [T^ (1971, P. 535) for the proof, and 
additional results regarding the convergence rate of the Edgeworth expansion): 


fni {x) ~ (x) 


1 - 
6 af ^ V 


Ci4 ! X 

+ I - 


24ar 


+ 




TTfi - 


X 


where the 6th Hermite polynomial is dehned as (x) = x® 
have /„. (x) ~ (x) (x) with: 

1 


72af 

+ 45x^ 


Cr,; 


(39) 
15. We hnally 


fn^ix) = 


y/^ai 


exp 
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= ( 1 + ^ 
“ ' 8af 
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5<3 \ 

24a,V 
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X -f 
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24cr? 
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V— 
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(- 


^2,4 
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V24af 24^1 V 
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—V-“ 

bi,4 


72d: 

bi,6 


(40) 


In this context, the likelihood maximization problem becomes: 
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exp 
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with 


^ ~ ' ^2 + :;2 + :;2 




(7f (7^ 


_ , hit h2t \ 2 

ht - 1 -y + -^ ) ^ 

U 1 U o 


'1 *^2 

Focussing for simplicity of exposure on symmetric noise distribution functions (therefore 
such that Cj^s = 0), we have: 


9ni (x) 9n2 (a^) = (^1,0 + bi,2x‘^ + &i,4a^^) (^2,0 + b2,2x‘^ + &2,4a:'^) 
= /3o + /32 x'^ + /34a:^ + + Ysx^ 


with; 


Yo = ^1,0^2,0 

1^2 = & 1 , 0 & 2,2 + & 1 , 2 ^ 2,0 

Ya = & 1 , 0 ^ 2,4 + ^ 1 , 2 ^ 2,2 + & 2 , 0 ^ 1,4 

/^e = &1,4&2,2 + &1,2&2,4 

/^8 = &1,4&2,4 

So we need to compute the following integrals, which can be obtained from the first 


27 



moments of the Gaussian distribution: 
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Finally, we have that: 
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(42) 

(43) 


We note that when ci ^4 = 02,4 = 0, that is when the third and fourth-order cumulant 
vanish, as would be the case for a Gaussian distribution, then we have A = 1, A = A = 
A = A = At = 0 , and we recover the maximum likelihood statistic of the Gaussian case: 


aG TT 

= max-exp 

a,a,,<72 Ct (7 i(T2 

In general, the presence of the additional terms implies a correction with respect to the 
Gaussian case. In the Gaussian case, one obtains the following explicit expressions for the 
variables involved in the maximization of the likelihood detection statistic H: 

2 ^2 _ /—2\ + 
a = a = [a ) 

2 ^2 _ /—2 ^ 2 \ + 

) 

Here, the expression for the likelihood detection statistic is more involved and since it is 
not clear whether any analytical solutions can be obtained for the values for a, ai, (J 2 , Ci^ 4 , 02,4 
that would lead to the maximum in Aml, one would need to resort to numerical optimization 
procedures. Taking the log, we have that: 
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C. Non-Gaussian Signal and Noise Distributions 


The methodology can also be extended to acconnt for the presence of deviations from 
the Gaussian assumption for both the signal and noise distributions. To do so, we use 
again the Edgeworth expansion to approximate the unknown noise distribution as /* {x) ~ 
f? (a^) 9 s (a:) and fn, (x) ~ f^. {x) (x) with: 
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In this context, the likelihood maximization problem becomes: 
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Focussing for simplicity of exposure on symmetric noise distribution functions, for which 
we have C3, ci,3, 02,3 = 0 , we obtain: 

gs (x) gni (x) gn2 (x) = (bo + b2X^ + b4X^) ( 61,0 + &l, 2 a^^ + &l,4a(^) (& 2,0 + b2,2X^ + &2,4a^^) 


= 7o + 72X^ + 743^^ + 76X^ + 78X^ + 783^^° + 78X^^ 
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with straightforward expressions for the 7 ^ terms as a function of the 62,1 and bi coeffi¬ 
cients. 

So we need to compute the following integrals, which can be obtained from the hrst 
moments of the Gaussian distribution: 
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Again, we need to compute higher-order moments of the Gaussian distribution, which 
is given by the following formula for a normally distributed variable X with mean p and 
variance: 


j=0 \ dj 
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where n\\ denotes the double factorial operator n\\ = ](([ (n — 2t) = n (n — 2) (n — 4)... with 

i=0 

For example, we have that: 
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Finally, we have that: 
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X (-^0 + At + -^21 + -^31 + I At + let + ht + hot + h2t) 


( 52 ) 


We note that when Cj ^3 = Cj ^4 = 0, that is when the third and fourth-order cumulant 
vanish for the noise distribution, we then recover the maximum likelihood statistic from 

ffl. 

V. CONCLUSIONS AND EXTENSIONS 

This paper analyzes the comparative efficiency of the standard CC detection statistic ver¬ 
sus an alternative detection statistic obtained by cross-correlating squared measurements in 
situations involving non-Gaussian noise (and signal) distributions and heterogeneous detec¬ 
tor sensitivities. We hnd that differences in detector sensitivities have a large impact on the 
efficiency of the CC detection statistic, which is dominated by the alternative statistic when 
these differences reach one order of magnitude. This effect is smaller in case of fat-tailed 
noise distributions, but it is magnihed in case noise distributions have skewness parameters 
of opposite signs. On the other hand, higher-order cumulants of the signal distribution do 
not have a material impact on the relative efficiency of the two detection statistics in real¬ 
istic situations where the signal is expected to be small compared to the noise. Since our 
methodology requires the estimation of higher-order moments/cumulants of the noise dis¬ 
tribution, we extend the maximum likelihood estimator to the case of non-Gaussian signal 
and noise distributions and manage to recover analytical expressions for the log-likelihood 
function in case these distributions can be approximated by Edgeworth-type expansions. 

Our methodology can be extended in a number of directions. We may first consider a 
setting involving a correlated noise component, typically regarded as environmental noise, 
in addition to the specihc instrumental noise. On a different note, we have considered so far 
colocated and coincident detectors, an assumption which would hold in the case of Einstein 
Telescope. On the other hand, our framework should be extended to apply to a network 
of separated detectors such as Advanced LIGO-Virgo detectors, or joint observations by 
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Advanced LIGO and Einstein Telescope. This extension is important because these are 
precisely the types of situations where differences in sensitivities are expected to be most 
substantial. 
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